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1. Introduction 



In recent years much effort has been devoted to the study of the inequivalence of en- 
sembles of systems with long-range interactions [H 121 El IH El E]- Systems where the 
interaction decays as ~ l/r"", with a less than or equal to the dimension of the system, 
are said to be of long-range. Such systems are non-additive, meaning that interface 
energies between macroscopic systems are not negligible compared to the bulk ener- 
gies. Non-additive, long-range interacting systems are found in many areas of physics, 
prominent examples being self-gravitating systems [3 E] and spin systems of mean field 
nature [21 IH E]- Some characteristics of long-range interactions have also been observed 
in driven systems with only local stochastic dynamics [10]. An important feature of sys- 
tems with long-range interactions is that thermodynamical properties can be dependent 
on from which ensemble they were derived [5]. For example, the specific heat in the 
canonical ensemble is always positive, but may become negative in the microcanonical 
ensemble [HI HI]. For more features, intricacies and examples of long-range systems we 
refer the reader to [3, [12] . 

It has been conjectured in [2] that a necessary condition for ensemble inequivalence 
to occur in statistical long-range interacting systems is that the system undergoes a 
first-order phase transition. This statement has not yet been proven generally, but has 
been put on a solid basis for many cases in [131 E] . 

A paradigmatic system with long-range interactions is the infinite-range p-spin model, 
also known as the Ising model with infinite-range p-body interactions. Notable variants 
are the Blume-Capel [15] and the Blume- Emery-Griffiths [16J models in the infinite-range 
limit , in which p = 2 and a spin can take the values Si = 1, 0, —1. It was shown in ^ 
and [1] that the Blume-Capel model and a model with infinite-range and short-range 
interactions exhibit ensemble inequivalence in the region of their phase diagram where 
first-order transitions occur. Also, a generalized XY model with two-body {p = 2) and 
four-body {p = 4) was studied in the microcanonical ensemble in [I7j. 

It has not been studied, however, what would happen if a long-range interacting system 
has randomness. This problem is of interest, because it is the target of extensive studies 
in the field of spin glasses [TBI [191 EO]. Almost all calculations so far have been done 
in the canonical ensemble for spin glasses with many-body interactions. It should be 
necessary to check if any different results may be derived in the microcanonical ensemble. 

It is quite difficult to analyze spin glass systems with many-body interactions using 
the microcanonical ensemble. Even the analysis in the canonical ensemble is plagued 
by the complications of replica symmetry breaking [TS], and we will have to develop 
ingenious methods to treat systems with randomness in the context of the microcanonical 
ensemble. We therefore set our goal more modest in the present paper and study the 
random-field Ising model with p-body interactions in the hope that some general features 
of the interplay between many-body interactions, randomness and first-order transitions 
may manifest themselves, leading to non-trivial consequences in ensemble inequivalence 
already in this relatively simple random system. 

This paper is organized as follows. In section |2] we introduce the model and derive the 
thermodynamic potentials, the free energy in the canonical ensemble and the entropy in 
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the microcanonical ensemble, and discuss some of their general features. In section |3] we 
compare the phase diagrams of both ensembles and discuss their differences. Our study 
is concluded in section HI 



2. The Model and General Considerations 

In this section we introduce the ferromagnetic p-spin model (= Ising model with p-body 
interactions) in random fields, derive the thermodynamical potentials in the canonical 
and microcanonical ensembles for later comparison in section [3] 

The ferromagnetic p-spin model in external fields describes a system of spins, in 
our case Ising spins {Si = ±1), with infinite-range p-body interactions and each site is 
subject to a random external field. It is defined by the Hamiltonian, 



with J > 0. Throughout this study we will set J = 1 without loss of generality. The 
external field is assumed to obey a quenched bimodal distribution given by 

P{h,) = ^6{hi~ho) + ^6{h, + ho). (2) 
2.1. Canonical Solution - Free Energy 

The partition function can be calculated in the standard manner by the introduction 
two auxiliary fields [18] 

Z = Tt j dmS (^Nm - ^ exp (^^NmP + /3 ^ h,S, 

= Tr y dm dm exp ^(3Nm^ — Nmffi + ffi Si + (3 hiSi^ 

= / drh dm exp {^l3mP - Nmm + ^ ln(2 cosh(m + I3h,^)^ , (3) 

where /3 is the inverse temperature. Evaluating this integral by the saddle point method, 
together with the saddle-point equations 



m = I3pm^ (4) 

N 

i ' 

m 



1 ^ 

— ^ tanh(m + (3hi), (5) 
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leads to the mean-field free energy per spin, 

I3f{m,l3) = /3(p-l)mP- ^^ln2cosh/3(pmP-^ + /ii) 

i 

= /3{p-l)mP -^\n2 cosh f3{pmP'^ + ho) 

-^ln2cosh/3(pmP-i - /lo), (6) 
where m obeys the self-consistent equation 

m = - tanh (3 {pm^'^ + ho) + - tanh (3 {pm^^^ — ho). (7) 

We have used the self-averaging property |18] (meaning that J^idi^i)/-^ equal to 
the configurational average of g{hi) for sufficiently large A^, where g is some physical 
quantity) to rewrite the site average in the first line of ^ by the randomness average in 
the second line. The free energy ^ shows standard behavior of first- and second-order 
transitions depending on the values of p and ho. Details of the canonical phase diagram 
will be elucidated in the next section. 



2.2. Microcanonical Solution - Entropy 




Figure 1: Microcanonical entropy per spin versus magnetization per spin, scaled to Sq, 
its value at m = 0, for p = 2, at external field strengths a) ho = 0.55 and 
b) ho = 0.95. In b) it is seen that the entropy for this case with a first-order 
transition is not a continuous function of the magnetization for energies lower 
than e', which lies below the critical energy for the phase transition, e' < Cc. 
The two maxima at e = ec, indicated by arrows, have equal height. 

The microcanonical entropy per spin is calculated from s = Infl/N, where fl is the 
sum of states for given values of m (magnetization per spin) and e (energy per spin). 
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The calculations are straightforward, as sketched in [A], and the result is 

, , , /l + m e + mP\ , 
s(m,e,hn) = - — In 
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In the microcanonical ensemble, the basic strategy, corresponding to the free-energy 
minimization in the canonical ensemble, is the maximization of the entropy. See figure [T] 
for the entropy as a function of the magnetization at different energies for p = 2. We 
have a paramagnetic phase, when the entropy at m = is the global maximum. When 
the transition is second order as in figure [l] a), the entropy flattens around m = as 
the energy decreases and develops a maximum at a non-zero value of the magnetization, 
while the region around m = becomes a minimum, which implies an instability of the 
paramagnetic phase. In the case of a first-order transition (figure [l] b)), as the energy 
decreases, the entropy develops a second maximum, and at a critical value ec the two 
maxima are of equal height. As the energy decreases further, this second maximum 
dominates. At some energy e' the entropy becomes a discontinuous function of the 
magnetization and the ergodicity is broken: Intermediate values of the magnetization are 
not allowed, and consequently, the system in the low-m region cannot evolve dynamically 
to the high-m region. As a consequence, the metastable state with m = has a divergent 
relaxation time ^4j. The energy, where ergodicity is broken, lies below the critical energy 
e' < ec for p = 2, but for p > 3 we have e' > ec for high ho, while for low ho the 
ergodicity breaks after the phase transition takes place. Figure [2] shows this property as 
the restriction on the accessible region in the (m, e)-plane, the non-convexity of the region 
being closely related to the ergodicity-breaking. This unusual behavior is not unique to 
this model, but has been observed in other systems with ensemble inequivalence [2llH|5]. 



3. Canonical and Microcanonical Phase Diagrams 

In this section we investigate the thermodynamic properties of the model (1) introduced 
in the previous section. We compare and discuss the phase diagrams for the cases p = 2 
and p > 2 in the canonical and microcanonical ensembles and show explicitly that in 
certain intervals of ho and e the specific heat of the ferromagnetic phase is negative in 
the microcanonical ensemble, which is a clear indication of ensemble inequivalence. 



3.1. p = 2 

The phase diagram of the model with p = 2, as shown in figure [3] a), is qualitatively the 
same as the phase diagram of the models discussed in [21 H] . For small ho there is a con- 
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Figure 2: The allowed magnetization per spin (shaded) versus the energy per spin for 
p = 3 at a) /iq = 0.5 and b) = 0.782. At low energy the allowed region 
of magnetization splits into two parts. The system in one region of magneti- 
zation cannot evolve to the other region under the microcanonical constraint 
of constant energy. The vertical dashed line shows at which energy ec the mi- 
crocanonical phase transition occurs. For ho = 0.5 the phase transition occurs 
before the ergodicity breaks, while for ho = 0.782 the ergodicity breaks before 
the phase transition. 



tinuous (second-order) transition in the canonical ensemble up to the canonical tricritical 
point, beyond which the transition becomes of first order [21]. In the microcanonical 
ensemble the continuous transition coincides, for small ho, with the canonical one, but 
persists to a higher value of ho up to the microcanonical tricritical point. The first-order 
transition in the microcanonical ensemble manifests itself as a jump in the temperature 
and is represented by two lines in the {ho, T) plane, the MC-m > and MC-m = lines 
in figure |3] a). The reason for this jump is that the entropy has two maxima, one at zero 
magnetization and one at finite magnetization indicated by the arrows in figure JT b). 
Since the temperature is calculated from the relation 1/T = {ds / de){e,m*) (see (19) in 



|Aj), and there are two values for m* , zero and another finite value, there is a temperature 
jump in the phase diagram due to {ds/de){e,m* = 0) 7^ {ds/de){e,m* 7^ 0) for ho > h^. 
Note that the upper microcanonical first-order line in the phase diagram in figure |3] cor- 
responds to m* 7^ (MC-m > 0) and the lower curve to m* = (MC-m = 0). Between 
those two lines, the phase of the system is not determined uniquely by the temperature 
and thus the two phases, ferromagnetic and paramagnetic, are not completely separated 
in the {ho,T) plane. 

Figure |4] shows the inverse temperature versus the energy per spin for several important 
values of ho in the region of the first-order transition. For ho > h^ = 0.8739, there is 
a small region where the specific heat is negative, shown as thick, red lines, around the 
transition energy. As ho grows, the temperature develops a jump at h^ ^ 0.8882, and at 
/iq = 0.908 the negative specific heat region disappears. Thus, the region in the {ho,T) 
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Figure 3: a) The canonical and microcanonical {ho,T) phase diagrams at p = 2. b) 
Details of the phase diagram, showing the region around the tricritical points, 
bounded by the box in a). The arrow points where the canonical spinodal line 
(C-S, red) touches the microcanonical-m > line (MC-m > 0, dash-dotted). 
For values of the external field strength < ho < H'q, the specific heat is 
negative for energies lower than the transition energy up to a certain energy 
value. At the temperature develops a jump at the transition energy. The 
red line touches the dash-dotted line as indicated by an arrow. 



plane where the system has negative specific heat is very narrow, < ho < h'o, as 
depicted in figure [3]b). 

A close look at figure [s] b) reveals that the canonical spinodal line (shown in red) 
touches the upper microcanonical line (shown as dash-dotted) [H] at ho = ^o(= 0.908), 
indicated by an arrow, where the specific heat becomes positive again. The dashed hori- 
zontal lines in figure |4] show the temperature where the canonical phase transition takes 
place. They can be obtained by the Maxwell construction |22] which was generalized to 
microcanonical systems in [25] . 

3.2. p = 3 

The microcanonical phase diagram changes drastically when many-body interactions 
come into play. 

Figure [s] a) shows the canonical phase diagram for p = 3. Here, nothing unexpected 
happens, the phases, ferromagnetic and paramagnetic, are cleanly separated by a first 
order transition. The microcanonical phase diagram as shown in figure |5] b) exhibits a 
far richer structure. The paramagnetic solution extends almost over the whole {ho,T) 
plane, except in a narrow region upper-bounded by the black line. Below that line, 
a ferromagnetic phase (FM) exists. However, the ferromagnetic solution extends up 
to the temperature indicated by red triangles. The blue dashed line is the canonical 
spinodal line. In regions M - I and M - II between the red triangles and the black line. 
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Figure 4: Inverse temperature versus energy per spin for p = 2. The thick, red hne 
indicates negative specific heat, because the temperature increases as the en- 
ergy decreases. The thick green dashed fines indicate tfie canonical transition 
temperatures, coinciding witfi tfie MaxweU construction. 



tfie ferromagnetic and paramagnetic solutions coexist, witfi tfie ferromagnetic solution 
fiaving positive specific fieat in region M - II, between tfie red and black lines. In tfie 
region M - I between tfie red line and red triangles, tfie ferromagnetic solution witfi fiigfi 
energy fias negative specific fieat wfiile tfie low energy solution retains positive specific 
fieat. At ho = H'q tfie region of negative specific fieat disappears and for Hq > h'^ all 
stable ferromagnetic solutions fiave positive specific fieat. 

To understand tfiis befiavior in more detail, we fiave to look at caloric curves T(e) 
at various values of tfie field strengtfi. In figure |6] we sfiow tfie temperature versus tfie 
energy for field strengtfis ho = 0.001, ho = 0.3, ho = 0.782 and ho = 0.999. At very low 
ho{= 0.001) tfie energy (12) is almost exclusively dominated by tfie magnetization and 
tfie paramagnetic solution exists only in a very narrow energy region (not sfiown). A 
notewortfiy feature of tfie ferromagnetic solution is tfiat at fiigfi energies it fias negative 
specific fieat (sfiown as tfiick red lines), wfiere tfie temperature increases as tfie energy 
is lowered. As tfie field strengtfi increases, tfie effect of randomness becomes more 
important and tfie ferromagnetic solution becomes destabilized in tfie fiigfi energy region. 
At ho = ho{= 0.782) tfie region wfiere tfie ferromagnetic solution fias negative specific 
fieat becomes metastable. At fiigfi field strengtfis, very close to one, almost tfie wfiole 
ferromagnetic solution is metastable in tfie sense described in tfie following section. 
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Figure 5: The canonical phase diagram a) shows standard behavior: the ferromagnetic 
(FM) and paramagnetic (PM) phases are separated by a first-order transition 
hne. The microcanonical phase diagram b) has a richer structure. Between 
the paramagnetic (PM) and the ferromagnetic (FM) phases he regions where 
both phases coexist. In M - I the ferromagnetic solutions with high energy 
have negative specific heat, while in M - II all ferromagnetic solutions have 
positive specific heat. The phase M - I disappears for > h'^. The dashed 
blue line is the canonical spinodal line. 



remaining only at the lowest energy values. At > 1, the randomness prevents any 
ferromagnetic order. 

3.3. Stability Analysis for j9 = 3 

To clarify what we mean by stable, unstable and metastable solutions, let us look at the 
entropy per spin as a function of the energy per spin. As stated earlier, the entropy has a 
maximum at zero magnetization. As the energy decreases, a second maximum emerges. 
These two branches of the entropy are in direct competition as shown in figure [T] a). 
Above the critical energy ec the maximum at m = is higher, and thus the paramagnetic 
branch wins, shown in black, while the ferromagnetic branch lies at lower entropy, shown 
red dashed. Below ec, the ferromagnetic branch wins (black), while the paramagnetic 
branch lies at lower entropy, shown red dash-dotted. Note that there is a sharp cusp at 
the transition, and the entropy is not smooth. This is a direct consequence of the fact 
that the unstable solutions (blue line) do not mediate the phase transition in bridging the 
gap between the metastable ferromagnetic and paramagnetic solutions. Thus, they can 
not smoothen the entropy around the transition point as is often observed in systems with 
long-range interactions including the van der Waals gas/liquid. The unstable solution 
corresponds to the minimum of the entropy with respect to the magnetization (see 
figure [S]). The caloric curves corresponding to the entropy branches in [T] a) are shown 
in figure [T] b). The paramagnetic solution (black) is stable above ec and metastable 
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Figure 6: At very low field strengths the energy is almost entirely determined by the 
magnetization, and the ferromagnetic solution dominates. At ho = 0.001 the 
paramagnetic solution is not shown. At high energies the specific heat is 
negative as indicated by a thick red line. As the field strength increases, 
the randomness destabilizes the ferromagnetic solution at high energies. The 
high energy branch for Hq = 0.3 and Hq = 0.782 represents the paramagnetic 
solution. At ho = /io(= 0.782) the whole of the negative-specific-heat solution 
is metastable. The ferromagnetic solution destabilizes completely at ho = 1. 



below (red dash dotted), while the situation for the ferromagnetic solution is the exact 
opposite: stable below (black) and metastable above (red dashed) ec- As the ergodicity is 
broken, the unstable solution (blue dashed) disappears. However, at high field strengths 
the ergodicity breaks before the phase transition takes place (figure [sjb)). Other models 
with many-body interactions, but without randomness, also show ergodicity breaking 
before the phase transition occurs [17]. 

The Maxwell construction cannot be used to bridge the gap between the microcanoni- 
cal and canonical ensembles. The conventional Maxwell construction yields the canonical 
transition temperature from the caloric curve /3(e), while the microcanonical Maxwell 
construction ^ yields the transition energy and the corresponding canonical transi- 
tion temperature. However, both procedures rely on the fact that the unstable solution 
interpolates the gap between the paramagnetic metastable and ferromagnetic metastable 
solutions. In our case, the Maxwell construction is impossible. For low ho, the unstable 
solutions do not bridge the metastable solutions. As shown in figure iTlb), the unstable 
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Figure 7: a) The different branches of the entropy for ho = 0.5. The black line is 
the equilibrium entropy, while the red dash-dotted line is the paramagnetic 
metastable solution and the red dashed line the ferromagnetic metastable solu- 
tion. The unstable solution is the blue dotted line. The green dotted, oblique 
line is the concave-cover construction, b) The corresponding caloric curves. 
Below ec the ferromagnetic solution is stable (black) and the paramagnetic 
solution is metastable (red dash-dotted), while above ec the paramagnetic 
solution is stable (black) and the ferromagnetic solution is metastable (red 
dashed). The unstable branch is shown blue dotted. 



solution (blue dotted) does not touch the paramagnetic metastable solution (red dash- 
dotted). On the other hand, for high Hq, the unstable solution disappears even before 
the phase transition occurs due to ergodicity breaking. For more details on the Maxwell 
construction we refer the reader to the literature [3 [221 ES] . 

It is useful to note here that the microcanonical entropy and the corresponding caloric 
curve in figure [T] can be reproduced from the generalized free energy (|6]). Details are 
delegated to [Bj 

Qualitatively similar results have been found for p > 3. 

3.4. The Limit j9 — > oo 

It is important to see what happens when we let p — )■ oo, because this limit plays 
important roles in some problems related to spin glasses [TSl [T^ I2U] . The canonical 
phase diagram can be obtained analytically. If we take the limit p — )■ oo in equations (|6| 
and ^ and assuming \m\ < 1, we see that m = is the only solution for |m| < 1 and 
the free energy for the paramagnetic phase is simply fp = — ln(2 cosh(/3/to))//3- In the 
same limit, p — t- oo, m = 1 is also a solution of ([T]) and the ferromagnetic free energy is 
fp = —1 irrespective of the temperature. Another solution with \m\ = 1, i.e. m = — 1, 
is not allowed for p odd as one sees easily from ([T]). For p even, m = —1 is also a solution 
but with the same free energy as m = 1. At the transition point the free energy should 
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Figure 8: A typical example of the entropy per spin versus the magnetization per spin 
for p = 3. a) Low field strength /iq = 0.5, at energies e = —0.25 (black), 
e = —0.3 (red), e = —0.34 (green dashed) and e = —0.4 (blue dash-dotted) b) 
High field strength = 0.782 at energies e = —0.49 (black), e = —0.55 (red) 
and e = —0.6 (green dashed). For ho = 0.782, the ergodicity is broken before 
the maximum at m > reaches the same height as the maximum at m = 0. 
The unstable solution corresponds to the minimum of the entropy. 



be continuous fp = fp, which leads to the condition for the critical line 

2 cosh l3ho = e^. (9) 



The canonical phase diagram is shown in figure [9j a) . 

In the microcanonical case, letting p — > oo in ([s]) and assuming \m\ < 1 leads to 
the conclusion that the only stable solution is at m = 0. With (19) the paramagnetic 
solution leads to 

e = —hotanh (3ho. (10) 

The ferromagnetic solution, \m\ = 1, yields e = — 1 and the temperature cannot be 
defined since the entropy has the unique value s = at the sole allowed energy e = — 1, 
which forbids us to evaluate ds/de. However, at zero temperature, the entropy of the 
paramagnetic solution has the same value, s = 0, as the ferromagnetic entropy. 

To better understand this behavior, we investigate the dependence of the microcanon- 



ical phases for large p in figure 10 a) p = 13 and h) p = 26. We see that the ferro- 
magnetic region becomes smaller as p increases, while the temperature region where the 
ferromagnetic phase coexists with the paramagnetic phase extends to higher and higher 
temperatures and the whole region in the phase diagram tends to be M - II. Conse- 
quently, the limit p oo gives different results than when the calculations are done 
with p = oo. 
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a) b) 




Figure 9: a) Canonical phase diagram for p = oo. The ferromagnetic (FM) and para- 
magnetic (PM) phases are clearly separated, b) Microcanonical entropy for 
p = oo at external fields ho = 0.99 (black) and ho = 0.5 (red). The system 
is paramagnetic at all temperatures, because e = (m = 0) has always the 
largest entropy. 
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Figure 10: Microcanonical phase diagram for a) p = 13 and h) p = 26. The regions where 
the pure ferromagnetic and paramagnetic phases exist shrink with growing 
p, while the mixed phases, M-I and M - II, blow up. 
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4. Summary and Discussion 



We have performed a comparison of the microcanonical and canonical phase diagrams of 
the ferromagnetic p-spin model with bimodally distributed random external fields and 
with infinite-range interactions. 

For p = 2, both ensembles coincide in the region of low values of the external-field 
strength ho, where the transition is of second order up to the canonical tricritical point. 
From there on the phase diagrams disagree. The microcanonical second-order tran- 
sition persists, off the canonical phase boundary, to a higher value of up to the 
microcanonical tricritical point. The first-order transition from the paramagnetic to the 
ferromagnetic state is accompanied by a temperature jump. The specific heat of the 
ferromagnetic phase, in the microcanonical ensemble, is negative in a very small inter- 
val of the external-field strength and only around the transition energy. All in all, the 
microcanonical phase diagram of the infinite-range model with random fields for p = 2 
is comparable to that of previously studied models [H H] . 

For the case p > 2, only a first-order transition exists in both ensembles. The two 
canonical phases are cleanly separated from each other, with the ferromagnetic phase 
lying, as expected, at lower temperatures than the paramagnetic phase. The micro- 
canonical case, however, gives completely different behavior. The paramagnetic phase 
extends down to low temperatures, while the ferromagnetic phase extends up to the 
canonical spinodal line below a certain value of the field strength Hq and up to some 
other temperature, below the canonical spinodal line. The mixed phase below Hq splits 
again into two regions, one where the ferromagnetic solution has positive specific heat 
and another region where the specific heat may take negative values. 

Also, no conventional Maxwell construction is possible for p > 2,m contrast to the case 
of p = 2, where it could reconcile the microcanonical and canonical phase diagrams. The 
reason for this is that, while for p = 2 the metastable solutions are joined smoothly by 
the unstable solutions, for p = 3 the unstable solution does not touch the paramagnetic 
metastable solution. For high values of the external field, the ergodicity is broken before 
the phase transition occurs, and the unstable solution does not even exist there. 

In the limit p ^ oo the canonical phase diagram can be obtained analytically. It is 
qualitatively comparable to the case of p = 3. In the microcanonical ensemble, however, 
only the paramagnetic phase exists at finite temperatures and the temperature of the 
ferromagnetic phase is ill-defined for p = oo, but for p large (i.e. in the limit p — )■ oo), 
all temperatures have both ferromagnetic and paramagnetic solutions. 

We have therefore found that subjecting the many-body Ising model to random ex- 
ternal fields yields many interesting and unexpected phenomena. Some already known 
features, like ergodicity breaking and negative specific heat, have shown up also in this 
model. New features of the system emerge when p > 2, such the impossibility of the 
conventional Maxwell construction or the breaking of ergodicity before or after the phase 
transition occurs, depending on the external control parameter. Clearly, these facts de- 
mand more detailed studies of the effect of randomness in long-range interacting systems. 
New microcanonical methods have to be developed to tackle the emerging difficulties 
that disorder brings to this system, especially in the context of spin glasses, and studies 
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in this direction are underway. 
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A. Microcanonical Entropy 



Here, we derive exphcitly the entropy, (|8j), from the sum of states and analyze its deriva- 
tives. 

Let us denote the number of sites with the field quenched to hi{= ho or —ho) and spin 
Si as Ns^hi- Then there are four possible denotations, since the field can be +ho or —ho 
and the spin can be ±1: A^++ for spins +1 at sites with hi = +ho, for spins —1 at 
sites with hi = +ho, and so on. Since half of the hi are +ho and the other half —ho, the 
number of states is given by 

(iV/2)! (iV/2)! 
A^++!iV_+ !A^+_!A^__!' ^ ^ 

The energy per spin is calculated from ([T| as 



+ ho- 



N " N 
/io(4r2++-m-l), (12) 



(13) 
(14) 
(15) 
(16) 





I-+ and used ^ 


Si/N 


n+ 


+ + n+. = ^(1 


+ m) 


n_ 


. + -- = ^(1 


— m) 


n+ 


1 




n+ 


1 

_ + n = -. 

2 





Using Stirling's formula together with (12) and the above relations, we can express the 
entropy per spin as a function of the energy and magnetization (|8|, 

/ 1 + m e + \ / 1 + m e + m*' 

^^"^'^'^"^ = -[—4 — i^J^H^ — ^ 

1 + m e + \ / 1 + m e + 
4 ^ Aho ) ^\ 4 ^ Aho 
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1 — m e + rri^ \ f \ — m e + mP 
In 



4 Aho J \ A Aho 

1 — m e + mP\^ f 1 — m e + mP\ , „ 

+ —7 — In ^ — + — -ln2. (17) 



4 Aho J \ A Ah^ 

The first derivative with respect to the magnetization is 

ds \ ( •pmP~^\ 1 / e + 

1 ; m - 1 + m 



dm 4 V /iq / 4 V /lo 





) In 


1 




4 




) In 


1 




4 



1 / pmP-^\ 1 / e 
— 1 + In - 1 + m + 



4 V ho ) A\ h 
' 1 + ^-- In- l-m + 



4 V ho J A\ ho 

1 / , pmP-^\ , 1 A e + mP 

— 1 m - 1 — m 



4 V /lo y 4 V ho 

for which m = is always a solution of ds/dm = 0. The first derivative with respect to 
the energy is the generalized temperature [S], 

,„) . -L i„ (1 - + + .n. - (. + 

^ ^ 4/io {l+m + {e + mP)/ho){l-m + {e + mP)/ho) 

from which the equilibrium temperature is obtained by setting m = m*, the value of the 
magnetization where the entropy has its global maximum. The second derivative of the 
entropy with respect to the magnetization at m = gives an indication of the order of 
the transition 



dw? 



p(l — p)mP '^f3{m, e) + 

1(1— pmP^^ /ho){m — (e + mP)/ho) 



+ 



2 1- {m- {e + mP)/ho] 



2 



l {l+pmP-'/ho){m + {e + mP)/ho) 

2 1 - (m + (e + mP)//io)2 ' ^ ^ 

Now, if we let m — )■ 0, 

d^s _ ~ 

- — -(e,m = 0)= lira p(l — p)mP ^/3(m,e). (21) 

am^ m-^o 

For j9 > 2 it is easy to see that [d"^ s / dim?){m — )■ 0+) < and the transition is always 
first order, because the paramagnetic state m = is always locally stable. 



B. Microcanonical Entropy from the Generalized Free 
Energy 

It is useful to notice that the microcanonical entropy can be obtained by taking the 
supremum of the temperature derivative of the generalized free energy /(m, /3), as defined 
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a) b) 

-n r-i 1 1 0.8 I \ 




-1 -0.8 -0.6 -0.4 -0.2 -1 -0.8 -0.6 -0.4 -0.2 



e £ 

Figure 11: a) The branches of energy as obtained from the generahzed free energy. Un- 
stable solutions are in blue, the low-energy branch is the ferromagnetic and 
the high energy branch is the paramagnetic solution, b) The microcanonical 
entropy derived from the generalized free energy is shown in black, while the 
green dotted line is the canonical entropy. The unstable branch is shown 
dashed. These branches correspond to the three possible solutions of the 
self-consistent equations ttlh. 



in equation ([6]), with respect to the magnetization. This is equivalent to the procedure 
discussed in [5]. 



s^,(e) = sup|/3^^(m*,/3)|, (22) 
where e = dPf{m*,l3) /d(3, is the energy obtained from the generalized free energy 



at the value of the magnetization determined by the above condition (22) and m* is 
determined by the solutions of equation ([t]). 

The canonical entropy is calculated from the canonical free energy 

/,an(/3)=mf{/(m*,/3)} (23) 

m* 

s...-^^% (24) 
where m* denotes all solutions of the self-consistent equations ([T]). In contrast to the 



calculation in equation (22), the temperature derivative is taken after the extremization 



with respect to the magnetization in equation (24). 



In figure 11 a) we show the three branches of the energy as the derivative of the 



generalized free energy for = 0.5. The paramagnetic and ferromagnetic branches are 
shown in black, while the unstable solutions are shown in blue. It is not indicated where 
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the ferromagnetic and paramagnetic branches are metastable. Comparing figure 
to figure IT] b) we note the following. When we calculate the canonical free energy, we 
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ignore the canonically unstable (blue in figure 11a)) solutions, as we look for the branch 
which minimizes the free energy. Such canonically unstable solutions correspond to the 
strictly convex part of the microcanonical entropy [24j. However, part of this canonically 
unstable solution is stable and metastable (ferromagnetic) in the microcanonical ensem- 
ble (black and red dashed in figure [T] b)). If, on the other hand, we keep track of all 
three branches of the generalized free energy, we can obtain the microcanonical entropy 



via relation (22) as shown in figure UTlb). Here, the black line is the microcanonical 



entropy, while the green dotted line shows the canonical entropy that would have been 



obtained from a conventional calculation as done in equation (24). 
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